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Abstract 

I briefly review the progress in perturbative QCD calculations in LEP era and 
discuss the prospects for the extension of their accuracy to higher orders. 

Higher order QCD calculations are essential in order to assess and precisely 
quantify our confldence in the Standard Model. Perturbative predictions at the 
Leading Order (LO) in the QCD coupling as rely on tree-level matrix elements 
and therefore provide only an approximate description of cross sections and distri- 
butions. The unphysical dependence on renormalization and factorization scales 
turns out to be quite large. 

The simplest higher order calculations are those for fully inclusive observables. 
Examples of these quantities are a{e'^e~ — > hadrons), the hadronic branching 
ratio of the Z and of the r lepton. For these quantities infrared (IR) singularities 
essentially cancel at the integrand level and thus accurate predictions exist up 
to next-to-next-to-leading order (NNLO) Moreover these observables are 
affected by small non perturbative corrections and thus they are particularly 
suitable for as measurements [^. 

The most detailed QCD studies at e+e^ colliders are based on event shapes 
and jet cross sections. Event shapes variables are quantities that characterize the 
structure of the hadronic event. Thus, with respect to fully inclusive observables, 
event shape distributions are useful not only to measure as but also to perform 
stringent QCD tests. The price to pay is that non perturbative effects are more 
important Q and that the perturbative calculation is more difficult. Real and vir- 
tual contributions have a different number of final state partons and they should 
be integrated separately in order to cancel IR singularities. This is only possible 
by combining analytical techniques with numerical (Monte Carlo) integration. In 
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practice one performs analytically only the part of the integration that produces 
the singularities. There are two methods of this kind: the subtraction method 
[Q] and the slicing method (5|0. The NLO calculation for all the relevant 3-jet 
observables was performed by using the subtraction method in Ref.0. 

At the beginning of LEP era the standard procedure was the comparison of 
the data with NLO QCD. Far away from the two-jet region the NLO calculation 
with the renormalization scale /U of order of the center-of-mass energy Q usually 
gives good fits to the data. On the contrary, very small (and unphysical) values 
of /u Q are required in order to extend the range of the fit. 

This is a signal of the fact that close to the boundary of the phase space (two- 
jet region) large logarithmic corrections of the form as log^ 1/y appear that spoil 
the perturbative expansion (here y denotes the variable that becomes small in the 
two-jet region). These contributions, whose origin is due to the emission of soft 
and collinear gluons, have to be resummed to all orders. The resummation for- 
malism was developed up to next-to-leading-logarithmic (NLL) accuracy for those 
shape variables that fulfill exponentiation [Q. This property is a consequence of 
two basic conditions: i) matrix element factorization] ii) phase space factoriza- 
tion. The first is a dynamical condition that relies on the general properties of 
soft and collinear emission. The second is a kinematical condition that depends 
on the particular quantity we consider. Both have to be satisfied in order to be 
able to perform resummation. The NLL resummed calculation is matched with 
the NLO prediction (NLL-hO(as)) to give a result that is everywhere as good as 
the NLO result and much better in the y ^ region. Extensive studies based 
on NLL-I-O(ag) calculations have been carried out over the past years ||2|. These 
studies have shown that resummed predictions are less sensitive to the renormal- 
ization scale /i and in particular enable us to avoid the choice of extremely small 
(unphysical) ^ in the two-jet region. 

A jet is qualitatively defined as a collimated spray of energetic hadrons, and 
is considered as a typical signal of parton dynamics at short distances. However, 
in order to perform quantitative studies, one needs a precise definition of jet. 
Once a jet definition has been chosen, the n— jet cross section is a function that 
depends only on the resolution parameter ycut- When ycut becomes small jet 
cross sections develop large as log^ 1 / ycut corrections similar to the ones that 
affect shape distributions in the two-jet region. As it happens for shape variables, 
the resummation of these corrections is possible if conditions i) and ii) above are 
satisfied. The old JADE algorithm |Q does not allow to perform the resummation 
since it induces strong kinematic correlations that spoil ii). The kx algorithm 
was introduced to make the resummation in the ycut region possible p!o| ]. 
NLL-I-O(ag) predictions exist for jet rates [0] and jet multiplicities [pTl . 

As far as NLO calculations are concerned, there were important developments 
in the last few years. In the first applications both the subtraction and the slicing 
methods were used with extensive partial fractioning on the QCD matrix element. 
This procedure is intrinsically process dependent. Later it was understood that 
both methods can be generalized in a process independent manner. The key 
observation is that the IR singularities can be singled out in a universal manner 
by using the factorization properties of soft and collinear emission. Today we 
have general algorithms that in principle allow to compute any observable at 
NLO, both with the subtraction |12, 13 1 and with the slicing method [0]. As 



complete numerical method was proposed in Ref. [||. 



far as LEP physics is concerned new calculations appeared for: 3-jet observables 
[13|], 3-jet with mass effects included [|l^] and 4-jet [p^]. 

With this situation one may wonder why we should do an effort to go to 
NNLO. There are several reasons to do that. The first one is that since reli- 
able perturbation theory starts at NLO, error estimate should start at NNLO. 
As far as LEP is concerned, the NNLO calculation of 3-jet observables would 
considerably reduce the error on as ||2|. Then we should consider that improved 
experimental techniques and higher luminosities will require a better control on 
the QCD background. 

A NNLO calculation requires in general three basic ingredients: i) NNLO 
parton distributions; ii) two-loop amplitudes; iii) knowledge of the IR behaviour 
of tree-level and one- loop amplitudes at ©(ag). 

In the last years enormous progress has been achieved in all these aspects. 
For a consistent evaluation of a NNLO cross section at hadron colliders NNLO 
(three loop) parton distributions are required. Even though their NNLO evo- 
lution kernels are not fully available, some of their Mellin moments have been 



computed and, from these, approximated kernels have been constructed 
Recently, the new MRST [^] set of distributions became available, including the 
(approximated) NNLO densities, which allows an evaluation of the hadronic cross 
section to (almost full) NNLO accuracy. 

Until two years ago no two-loop amplitude depending on more than one scale 
was available. Thanks to the calculations of important two-loop master diagrams 



[20 1 and to the great progress in the reduction of tensor integrals the first 



calculations of two-loop 2—5-2 amplitudes have recently appeared in QED 



and in QCD |23|. The IR singularities appearing in these amplitudes are in 



agreement with the general prediction in Ref.|24|. These results will be relevant 
for the calculation of Bhabha scattering and two-jet cross sections in hadron 
colhsions at NNLO. 

To perform a NNLO calculation one has to combine the two-loop amplitude 
with the one-loop correction where one parton is unresolved and with the tree- 
level contribution where two partons are unresolved. The kernels that control 



soft ^] and collinear ||2^, 27 1 singularities appearing in one- loop amplitudes 
have been computed. The IR singularities appearing in tree-level amplitudes are 
more complicated because many soft/collinear limits have to be considered. All 
these limits have been studied pSl, pOl and the corresponding kernels have been 



computed pl-P] 



The step that remains to be performed is to combine all these ingredients 
to construct general algorithms to handle and cancel infrared singularities. This 
step is more difficult than at NLO since the pattern of IR singularities is much 
more complicated. Nevertheless some applications where some progress in this 



direction has been achieved recently appeared |33|, |34|. 



The results of Refs.[25|-|31] are relevant not only to perform NNLO calcula- 
tions, but also to extend the accuracy of resummed calculations at NNLL. Trans- 
verse momentum (kx) distributions of high-mass systems (lepton pairs, vector 
boson, Higgs...) in hadronic collisions are affected in the small kT limit by large 
logarithmic contributions of the same (infrared) nature of the ones present in 
event-shape distributions in the two-jet limit. 
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In Ref ■ [^] the structure of these large corrections was studied at 0{a'^) up 
to NNLL accuracy. This calculation was performed in a general (process inde- 
pendent) manner by exploiting the universal nature of these corrections. The 



results of Refs.[25|-p(| were used to construct improved approximations of the 
relevant matrix elements that allow to control all the infrared singular regions 
responsible for the appearance of the logarithmic contributions. This method, 
even if strongly dependent on the special kinematics of this class of processes, 
could be extended in the spirit of Ref.[13| to more general cases f\. 

In Ref. 1 34] the calculation of the soft and virtual corrections to Higgs boson 
production at hadron colliders was presented. This calculation was done by com- 
bining the recent results |35] for the two- loop amplitude gg ^ H in the large 
mtop limit with the soft factorization formulae for tree- level |^, 29 1 and one- loop 
[25, 26| amplitudes 0. From the theoretical side this calculation is very important 
since it provides a check of the cancellation of the IR poles from to 1/e be- 
tween real and virtual contributions. From the phenomenological side the results 
give a first consistent estimate of the QCD corrections to this important process 
at NNLO. 

Up to a few years ago NNLO calculations, if doable, were considered very far 
in the future. With the progress achieved in the recent years we can be more 
confident that these calculations will be feasible in the LHC era. 
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